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Abstract. Fermi coordinates are the natural generalization of inertial Cartesian 
coordinates to accelerated systems and gravitational fields. We study the motion of 
ultrarelativistic particles and light rays in Fermi coordinates and investigate inertial 
and tidal effects beyond the critical speed c/a/2- In particular, we discuss the black- 
hole tidal acceleration mechanism for ultrarelativistic particles in connection with a 
possible origin for high-energy cosmic rays. 
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1. Introduction 



Physically meaningful interpretation of the measurement of relative motion from the 
point of view of an accelerated observer in a gravitational field requires the introduction 
of a special coordinate system (i.e. Fermi coordinates) along the worldline of the 
observer. In this coordinate system, the equations of relative motion reveal tidal 
and inertial effects for ultrarelativistic motion (with speed exceeding the critical value 
c/ \/2) that are contrary to Newtonian expectations. While there are general treatments 
of inertial and tidal effects in Fermi coordinates ^1 |21 El E] and the special case of 
ultrarelativistic motion of particles has been discussed in our recent papers El EI, 
the purpose of this work is to present a more systematic and complete description of 
motion beyond the critical speed c/ \/2. To this end, the equations of motion in Fermi 
coordinates are discussed in section El and inertial effects are considered in sections El IH 
and appendix A. Section El is devoted to tidal effects. A summary and brief discussion 
of our results is contained in section El 



2. Equations of motion in Fermi coordinates 

Imagine an accelerated observer in a general relativistic spacetime following a worldline 
x^(r), where r is the proper time along its trajectory. The local axes of the observer 
are given by an orthonormal tetrad frame ^^(^a) ^^^^ carried along its path according 
to 

where $(a)(/3) is the antisymmetric acceleration tensor, A'^^g-j = dx^/dr is the local 
temporal direction and A'^^j-j, i = 1,2,3, form the local spatial triad, and where (here 
and throughout this paper) Greek indices run from to 3, Latin indices run from 1 to 

3, the signature of the metric is +2 and c = 1, unless specified otherwise. 

In analogy with the electromagnetic Faraday tensor, we decompose the acceleration 
tensor into its "electric" and "magnetic" components 4>{a)(p) (~ci, c*;); these are given 
by the translational acceleration a(r) and the rotational frequency tL>(r), respectively. 
That is, the acceleration of the reference trajectory is given by a*A^(-.-) and a; is the 
frequency of rotation of the spatial triad with respect to a local nonrotating (i.e. Fermi- 
Walker transported) triad. 

Let us next establish a Fermi coordinate system (T, X) in a neighborhood of the 
reference worldline It turns out that Fermi coordinates can be assigned uniquely 
only to spacetime events that are within a cylindrical region of finite radius along 
the observer's worldline. At each event x^{t) on the observer's worldline, consider 
all spacelike geodesic curves that are normal to this reference worldline at this event. 
Each event x'^ on the resulting hypersurface and within the cylindrical region under 
consideration is connected to x^{t) by a unique spacelike geodesic curve that starts at 
x^(r) and whose tangent vector ^'^ at this event is normal to the reference worldline 
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(that is, ^/iA^^Q-j = 0). The event is assigned the Fermi coordinates = (T, X), 
where T = t and 

X^ = aeA/). (2) 

Here a is the proper length of the normal spacelike geodesic segment connecting x^{t) to 
x'^. The Fermi coordinate system is admissible in a cylindrical spacetime region around 
x^(r) with radius |X| ~ 7^, where TZ is the minimum radius of curvature of spacetime 
along the reference worldline |2]. 

The observer, at each instant r = T of its proper time "sees" a locally Euclidean 
three-dimensional space and "instantaneously" determines distances within it using the 
(spacelike) geodesic lengths beginning from its position at the spatial origin of the 
new Fermi coordinates. A nearby particle worldline that punctures the sequence of 
three-dimensional spaces at X = (X(T), Y{T), Z{T)) is thus a graph over the reference 
worldline and each point on the graph is connected to the observer by a "locally straight" 
line (i.e. a geodesic) that is normal to the observer's worldline. From the viewpoint of 
the observer, the particle has relative coordinate velocity V = dX/dT and relative 
coordinate acceleration dV/dT. 

The general equation of motion for a test particle of mass m in the Fermi frame is 

+ r^«-T— T- = ^^ (3) 



ds2 ds ds 

where s is the proper time along its worldline, — ds^ = ^'^^.(T, X)dX^dX'', and 
mA^ = is the external force acting on the particle. 
Equation Q can be written as the system 

d^T . dX°dX^ . 
d^X' ■ dX"dX^ _ . 
Using the identity 

d^^d^^^ d^ (6) 

ds2 ^ dT2 ' ^ ' 

and the Lorentz factor of the particle 

r(r.x):=^, (7) 

equation (0) can be expressed as 

!i^ + (ri,-r»,n^^^^ = i(^'-AV'). (8) 

Let := dX^/ds. We note that = r(l, V) and the physical content of the 
equation of motion Q is contained in equation (jS)) together with Ufj,U^ = —1 and 
Ufj^A^ = 0. These can be written respectively as 

1 



^-goo - 2goiV' - gijV'V^ 



(9) 
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and 

^00 + goiV' 

It is crucial to recognize that in Fermi coordinates the velocity V satisfies the condition 
I V| < 1 only at X = 0; indeed, away from the reference trajectory | V| could in principle 
exceed unity in accordance with the equation of motion (|H)). 

3. Inert ial effects 

It is natural to begin our discussion with an accelerated observer in Minkowski 
spacetime. In this case, spacelike geodesic segments are straight lines; therefore, 
the construction of Fermi coordinates is simple. In fact, equation Q reduces to 
— x^{t) = X*A'^^j^(r). Differentiating this relation we find that 

dx^ = [PA^(o) + Q'^\j)]^X' + X'^,^dX\ (11) 
where P and Q are given by 

P(T,X) = l + a(T) -X, Q(T,X) =a;(T) X X. (12) 
The components Qf^^, of the Minkowski metric tensor in Fermi coordinates are given by 

goo = -P^ + Q^ Qoi = Qi, Qij = Sij. (13) 
Moreover, det((7^j/) = — and the inverse metric is given by 

^" = -^ 9'' = % 9- = k-^- (14) 
The Christoffel symbols can be evaluated using ()13|) and ()14|) : the nonzero components 



are 



J- 00 - p ' ■'-oj-p' J- oi - -l^ufctu +— ^j, [lb) 



rj)o = Pai - ^^Qi + [a; X (a; X X) + X X],, (16) 

where an overdot denotes differentiation with respect to T and S is defined by 

S(T) := a + ax a;. (17) 

In Minkowski spacetime, the Fermi coordinates are admissible for P^ > Q^. We 
note that the structure of the boundary region P^ = has been discussed in detail jH]. 
In general, the boundary of the admissible region at a given time T is a real quadric 
cone. If a; = 0, this surface degenerates into coincident planes. Moreover, if a; 7^ 0, 
but a • a; = 0, the boundary surface is a hyperbolic cylinder for Icjp < |ap, a parabolic 
cylinder for \uj\'^ = |ap, and an elliptic cylinder for |a;p > |ap. As is well known, for 
a = 0, the boundary surface is a circular cylinder of radius 
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We now consider the equations of motion in Fermi coordinates. Equations Q 
and (fTUj) are given by 

l = p2_(Q + V)2>0, (18) 
^ (Q + V)-A 

Hence, equation (jH)) has the form 
d^X 

+ 2a;xV + a;x(cjxX) + cjxX 



dT2 

+Pa-475(Q + V)(S-X + 2a- V) =:F, (20) 
where J- represents the external force per unit mass 

:f-^\k + v] (21) 

•^"P^^ C2(P2_Q2_Q.V)^J- ^^^^ 

It is interesting to note that in equation piU) the purely rotational inertial accelerations 
are essentially the same as in the nonrelativistic theory [HlITj. 

Inertial accelerations have been discussed by a number of authors [SI IH HI CHI HH 
EIEJEIE]- In particular, it has been shown [7] that for a; = and a = 0, the inertial 
acceleration experienced by the particle parallel to its motion is given to lowest order in 
a by — a - V(l — 2V^/c^), where V = V/V is the unit vector tangent to the spatial path 
of the particle; therefore, there is a sign reversal for V >'V^ = c/\/2 with consequences 
that are contrary to Newtonian expectations. 

Note that the acceleration A and the inertial acceleration a occur in rather different 
ways in equation ()20|). This difference is a consequence of the absolute character of 
acceleration in the theory of relativity; that is, the fact that an observer is accelerated 
is independent of the choice of coordinates. In particular, while the critical speed 
associated with A is c (on the basis of equation (j2II)), the critical speed associated 
with a in equation is c/ \/2. 

To obtain the limiting case of lightlike motion, we let ds = m dA, where A is an affine 
parameter. In the limit as m — > 0, we have that = (Q + V)^; hence, equation (pUll 
is valid with ^ = (see appendix A). 



4. Rindler observer 



For a Rindler observer ^H] (that is, an observer in hyperbolic motion in Minkowski 
spacetime), consider an inertial frame in which the Rindler observer is at rest at 
(0, 0, 0, Zq) and has uniform translational acceleration g along the z-axis. The observer's 
worldline is given by 

t = -sinh5fr, x = ?/ = 0, z = zq-\ — (— 1 + cosh^fr). (22) 
9 9 
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The nonrotating orthonormal tetrad frame along this worldhne has the form 

A'^(o) = 7(1, 0, 0, (3), A^(^ = (0, 1, 0, 0), (23) 
A'^(2) = (0,0,1,0), A^(3)=7(/?, 0,0,1), (24) 

where f3 = tanh^fr and 7 = cosh^fr, and the transformation from Fermi to inertial 
coordinates is 

t= (^Z + ^^ sinh gT, x = X, y = Y, (25) 

z = - ^ + + ^) cosh gT. (26) 

The Rindler coordinates are admissible for T,X,Y G (—00,00) and Z G (—1/(7,00). 
We note that the manifold where Z = —1/g is the Rindler horizon. It is related 
to the null cone in the inertial frame, since equations ()25|) and (j^U]) imply that 
{Z + l/gf = {z-Zo + l/g)'-t\ 

Imagine a free particle moving along the 2;-axis with speed Vq, —1 < Vq < 1, 
according to 

z = vot + zq. (27) 
From the viewpoint of the Rindler observer, this motion is given by 

Z=-(-l+ ^ ^ ^ . ^ J . (28) 



g \ cosh gT — vq sinh gT^ 

For fo < 0, the particle descends monotonically toward the horizon. For > 0, its 
path crosses the path of the Rindler observer {Z = 0) at two events corresponding to 
T = t = and T = Ti, where 

T. = il„(i±i^). (29) 

This value corresponds to inertial time ti = 2'yQVQ/g, where 70 is the Lorentz factor 
associated with vq. At either event 
d^Z^ 



dT2 



= -g{l-2v',). (30) 
z=o 



On the other hand, Z{T = 0) = vq and Z{T = Ti) = —vq. That is, the free particle 
ascends from Z = to Zq = (70 — l)/(7 at Tq = Ti/2 such that tanh^fTo = Vq, then 
descends back to Z = at Ti and finally approaches the horizon Z = —1/g a.s T 00. 
For Vq = 1/ \/2, the particle has no acceleration at the crossing events. 

We note that for ultrarelativistic motion beyond the critical speed {vq > l/\/2), 
the inertial acceleration of the free particle has the opposite sign at crossing events 
compared to intuitive expectations based on Newtonian mechanics; but, other aspects 
of the motion are not affected. For instance, the critical speed does not appear to play 
a role in the proper temporal intervals between the crossing events. The relevant proper 
time of the Rindler observer Ti given by equation (^^1) is always less than the proper 
time of the free particle ^1/70 = 270^0, since for < fo < 1, 

In f 1±^) < ^2^. (31) 
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That is, the length of the geodesic segment is maximum. This circumstance, usually 
called "the twin paradox," is a reflection of the absolute character of acceleration in the 
theory of relativity. 

For the limiting case of a null ray {vq = ±1), P = 1 + gZ = ±Z = exp{i:gT) 
and Z = gexp{±gT). The ray thus ascends monotonically toward infinity or descends 
monotonically toward the horizon always with positive acceleration. Moreover, the 
Fermi coordinate speed of the ray can range from zero to infinity. 



5. Tidal effects 

We now turn our attention to motion in a gravitational field PJEllSlll]- The spacetime 
metric in Fermi coordinates then takes the form 

^700 = -P' + ^RoiOjX'X^ + 0(|Xn, (32) 

go^ = Q^-l ^RojikX'x^ + 0(|xn, (33) 



g,^ = 5,, - ^RikjiX^X' + 0(|Xn, (34) 



where 



^R(xp-,5{T) — Rfj,upa^^ (^a)"^" (13)'^'' i-y)'^'^ (S) (35) 

is the projection of the Riemann tensor on the tetrad frame of the observer. 

To simplify matters, we concentrate on tidal effects only and assume that the 
observer follows a geodesic along which a Fermi coordinate system is constructed based 
on a parallel-propagated tetrad frame. Moreover, we neglect external forces on the test 
particle. In this case, the equations of motion of the free particle in the Fermi system, 
i.e. equations (jH|)- (fTm) . reduce to 

^ + ^RoiOjX^ + 2^R,kjoV''X^ 

+ (^2^Rok,oV'V'' + '^^R^k.eV'V' + '^^RokjiV'V'^V'^ + 0(|xn = 0, (36) 



and 



l- = l-V^ + ^RmjX'X^ + ^^R^,,i,X^V'X^ 



+ ^-'^R^k,lV'X''V^X' + 0(|Xn > 0. (37) 

Equality holds for null rays (F = oo). In this case the right-hand side of equation (jHH) 
is a first integral of the differential equation ()36|) and higher-order tidal terms cannot be 
neglected. The Fermi coordinate system is admissible within a cylindrical region along 
the observer's worldline with |X| < IZ, where IZ^'^iT) is the supremum of \^RaP'ys{T)\. 
Equation (jHBj) is the geodesic deviation equation in Fermi coordinates, since it represents 
the relative motion of the free particle with respect to the fiducial observer that has been 
assumed here to follow a geodesic. 
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To illustrate the tidal effects for test particles and null rays, we consider the 
gravitational field of a Kerr black hole. The Kerr metric, in Boyer-Lindquist coordinates 
(t, r, {}, 0), is given by 

- ds^ = -de + dr^ + dd^) + (r^ + a') sin^dd<p' 

+ 2GM^{dt- asin^^dcpf, (38) 

where E = + cos^ i} and A = — 2GMr + a^. The Kerr source has mass M and 
angular momentum J = Ma. For a black hole, the specific angular momentum a is such 
that < a < GM. The observer is taken to be on an escape trajectory along the axis 
of rotation: It starts at r = with tq > y/Sa on this axis and moves radially outward 
reaching infinity with zero kinetic energy. Its geodesic worldline is given by 



dt _ r^ + a^ dr _ 2GMr 

d^~ r^- 2GMr + a? ' d^~ \l r^ + a^' ^ ' 

To determine Fermi coordinates, we choose a nonrotating orthonormal tetrad frame 
along the observer's worldline such that X^^^^ is parallel to the z-a.xis. In {t, r, -(9, 0) 
coordinates, we have 

A'^(o) = (t,f,0,0), A% = (ir, 1,0,0), (40) 

where t = dt/dr and r = dr/dr are given in display (jHOl)- Using the axial symmetry 
of the spacetime, A''^-^^ and A'^^2) can be chosen uniquely up to a rotation about the 
2;-axis. Once a pair is chosen at r = 0, then A'^^^-j and X'^^^) parallel propagated along 
the reference worldline. We do not require the explicit transformation from Boyer- 
Lindquist to Fermi coordinates to obtain the curvature components that appear in the 
approximate equations of motion (jH?)|) - ()37|l . Indeed, the symmetries of the Riemann 
tensor make it possible to represent these quantities as elements of a symmetric 6x6 
matrix with indices that range over the set {01, 02, 03, 23, 31, 12}. The result is 



E B 
B -E 



(41) 



where E and B are 3x3 matrices that are symmetric and traceless and correspond 
respectively to the electric and magnetic parts of the Riemann tensor. Moreover, 
E = diag(— fc/2, —k/2, k) and B = diag(— g/2, —q/2, q), where 

k = -2GM%^^, (42) 
3r2 - n"^ 

q = 2GMa (43) 
[r^ + a^y^ 

The curvature components are computed along the worldline of the observer, i.e. at 
(T, 0) in Fermi coordinates. The electric and magnetic curvatures in (I42p and (j43p . 
respectively, only depend upon r. Therefore, it suffices to integrate the equation for r 
in (j39|) to obtain r(r), which when substituted in (I42|) and (j43p with r ^ T results in 
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Figure 1. Graphs of Fji defined in 148|) versus T/{GM) approximated by numerical 
integration of equation (gTJ with initial Fn equal to 5, 10, 20, 50, 100 and 1000. We 
assume that at T = 0, ro = 20 GM and Z{Q) = 0. In these graphs a — GM; however, 
the corresponding graphs with a = are indistinguishable from those given here. 



k{T) < and q{T) > 0, respectively. The equation of motion ()3()|1 can thus be written 
as the system 



X - -kX (l-2X^ + -Y^ - -Z^ 
2 V 3 3 



+ -kX{5YY - 7ZZ) 



Y 



+ q[XYZX - ZY{1 + X^ 
1 



2YZ] = 0, 



-kY 1 



2Y^ + ^X^ - '^zA + ^kY{5XX - 7ZZ) 



- q[XYZY - ZX{\ + Y^) - 2XZ] = 0, 



Z +kZ 



2Z^ + -{X^ + Y^) 
3 



0, 



+ -kZ{XX + YY) 



(44) 



(45) 



(46) 



- q{XY - XY){1 - Z^] 

where we have neglected higher-order tidal terms for the sake of simplicity. These 
generalized Jacobi equations reduce to the standard Jacobi equations when the relative 
motion is so slow that the velocity terms can be neglected. 

Consider, for instance, the one-dimensional motion of the test particle along the 
^-direction. Equations with X(T) = Y{T) = for T > reduce to 

^ + Ml - 2Z')Z = 0. (47) 

The behavior of the solutions of this equation, which contains the critical speed 
K = 1/72, has been studied in detail in our previous work Suppose that 
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Figure 2. Graphs of Fj^ defined in H5()|l versus T/{GM) approximated by numerical 
integration of equation (|49|l with a = GM and initial rj_ equal to 5, 10, 20, 50 and 
100. We assume that at T = 0, tq = 20 GM and X(0) = 0. Note that T±_ diverges at 
T/{GM) w 23.7, 13.1, 6.5 and 4.0 for the initial Lorentz factors T ^{Q) = 10, 20, 50 
and 100, respectively 



Z{T = 0) = and Z{T = 0) > l/-\/2; then, the particle starts from the observer's 
location and decelerates along the z-axis toward the critical speed 1 / \/2. In this case F 
given by equation (jHTj) reduces to 

Til = , (48) 

Vl - Z2 + kZ^ 

whose behavior is depicted in figure H This extends our previous results [Hj and is 
consistent with the rapid decrease of curvature k away from the black hole. 

A more interesting situation arises for motion normal to the direction of rotation 
of the black hole. Equations (jlll)-(jini) with Y{T) = Z{T) = for T > reduce to 

^ - - 2X')X = 0, (49) 

which turns out to be the general equation for one-dimensional radial motion 
perpendicular to the rotation axis of the black hole as a consequence of the axial 
symmetry of the configuration under consideration here. The corresponding Lorentz 
factor is given by (jTTj) as 

1 



X2 - iA;X2 



(50) 



We note that for ultrarelativistic motion, the particles gain tidal energy and reach the 
speed of light as in figure |21 for an extreme Kerr black hole (a = GM) and in figure 121 for 
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Figure 3. Graphs of r_L defined in H5()|l versus T/{GM) approximated by numerical 
integration of equation H49(l with a = and initial V± equal to 5, 10, 20, 50 and 100. 
We assume that at T = 0, ro = 20 GM and X(0) = 0. Note that T ±_ diverges at 
T/{GM) w 23.5, 13.0, 6.4 and 3.9 for the initial Lorentz factors T ^{Q) = 10, 20, 50 
and 100, respectively. 

a Schwarzschild black hole (a = 0). It follows from a comparison of figures |21 and El that 
the tidal acceleration mechanism does not depend sensitively upon the specific angular 
momentum of the black hole; in fact, it appears that the time it takes for Fj^ to diverge 
is in general slightly longer for a = GM as compared to the a = case. Moreover, the 
divergence of T± indicates that our test-particle approach breaks down. 

In two recent papers fBl, we approximated r_|_ by (1 — X^)^^/^ and plotted its 
approach toward infinity as a consequence of the ultrarelativistic tidal acceleration 
mechanism. The present calculation of r± is based on an improved approximation 
scheme and goes beyond our previous work. To reach the speed of light — or, 
equivalently, for the timelike particle worldline to become null — the Fermi speed X 
should reach a value that is greater than unity as k in equation (jSU)) is negative. We 
expect that this value of X is reached at a later time T as demonstrated by a comparison 
of figures m and El with our previous results [Hj. Even the present calculations are 
approximate as they are based on dropping the higher-order tidal terms. It is clear 
that at a sufficiently long time T ~ 7^, the Fermi coordinate system breaks down. 
The kinematic breakdown of Fermi coordinates, i.e. their inadmissibility, is logically 
independent of the divergence of F^, which indicates the failure of the dynamical theory 
presented here. Our test-particle approximation ignores any back reaction; furthermore, 
for charged particles the electromagnetic field configuration near the black hole could 
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enhance or hinder the tidal acceleration mechanism. In any case, this breakdown of 
our dynamics nevertheless implies that sufficiently energetic particles emerge from the 
vicinity of the black hole after having experienced tidal acceleration by the gravitational 
field of the collapsed source. In this way, highly energetic particles may be created by 
microquasars in our galaxy. It is important to emphasize that the tidal acceleration 
mechanism is independent of the horizon structure of the black hole. These results 
are interesting in connection with the origin of the highest energy cosmic rays since 
cosmic rays with energies above the Greisen-Zatsepin-Kuzmin limit (~ 10^° eV) are not 
expected to reach the Earth from distant galaxies [T7] . 



6. Discussion 



Fermi coordinates constitute a geodesic coordinate system that is a natural extension 
of the standard Cartesian system and is indispensable for the theory of measurement 
in relativistic physics. We have discussed the general equation of motion of a pointlike 
test particle, as well as the limiting case of a null ray, in Fermi coordinates. Inertial and 
tidal effects of ultrarelativistic particles with speeds above the critical speed c/ \pl ~ 0.7c 
have been emphasized. This work goes beyond our previous work in this direction and 
strengthens the basis for our results and conclusions OIHIIZI- 



Appendix A. Null rays 

Let us compute the effective external force per unit mass given by equation for the 
electromagnetic case, namely, 

= ±F^,^ (A.l) 
m ds 

in accordance with the Lorentz force law for the motion of a particle of charge q in an 
external field F^^,^. Neglecting radiation reaction, we find that 

r = q (n + F',V^ - F'^^V^Vn, (A.2) 

where A is the affine parameter defined hj ds = mdX. 

A massless charged particle does not exist; therefore, we assume that g — > as 
m — s> 0. With this assumption, it turns out that the massless limit of the trajectory of 
a charged particle is a null geodesic. 

There exist nongeodesic null curves in Minkowski spacetime. These may be 
physically interpreted as follows: Imagine the path of a null electromagnetic ray that is 
reflected from a collection of mirrors fixed in different places in space. The path consists 
of straight (i.e. geodesic) null segments separated by mirrors. Next imagine a limiting 
situation involving an infinite number of such idealized pointlike mirrors. The resulting 
smooth curve is a nongeodesic null curve. An example of such a curve is 

x^{X) = (aA, a sin 6* sin A, a sin 6' cos A, aA cos ^^), (A. 3) 

where a and 6 are constants. The generalization to curved spacetime is straightforward. 
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